Uniaxial contribution to the magnetic anisotropy of La .67Sro.33Mn03 thin films 
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Lao.67Sr33Mn03 (LSMO) thin films under compressive strain have an orthorhombic symmetry 
with (110)o and (001) o in-plane orientations. (The subscript o denotes the orthorhombic symme- 
try.) Here, we grew LSMO on cubic (LaA103)o.3-(Sr2AlTaOe)o.7 (LSAT) substrates and observed 
a uniaxial contribution to the magnetic anisotropy which is related to the orthorhombic crystal 
structure. Since the lattice mismatch is equal in the two directions, the general understanding of 
anisotropy in LSMO, which relates the uniaxial anisotropy to differences in strain, cannot explain 
the results. These findings suggest that the oxygen octahedra rotations associated with the or- 
thorhombic structure, possibly resulting in different Mn-O-Mn bond angles and therefore a change 
in magnetic coupling between the [110] o and [001] o directions, determine the anisotropy. We expect 
these findings to lead to a better understanding of the microscopic origin of the magnetocrystalline 
anisotropy in LSMO. 
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I. INTRODUCTION 

The perovskite oxide Lai_ x A K Mii03 (A=Ca, Ba, 
Sr) has initiated a substantial body of research due 
to its colossal magnetoresistance 1,2 . Extensive theo- 
retical studies and experimental investigations utiliz- 
ing Lai-^A^MnOs perovskites in bulk form revealed 
a strong coupling between lattice distortions and mag- 
netism, which substantially modify magnetic proper- 
ties such as magnetoresistance and Curie temperature 3 -^. 
Lao.67Sr33Mn03 (LSMO) has the highest Curie tempera- 
ture (370K) and a 100% spin polarization^. LSMO can 
be coherently grown on a range of commercially available 
perovskite substrates, such as e.g. NdGaOs (NGO) and 
SrTi03 (STO). The epitaxy stabilizes a different crystal 
structure which modifies the magnetic properties. Espe- 
cially magnetic anisotropy is shown to be very sensitive 
to the LSMO crystal structure 7 - - — . When anisotropic 
strain is applied to the LSMO the magnetocrystalline 
anisotropy becomes strongly uniaxial—^, which is a use- 
ful tool to tailor the magnetic properties for device ap- 
plications. 

In the case of isotropic tensile strain, e.g. tetrago- 
nal LSMO thin films on cubic STO (001) c substrates, 
the magnetocrystalline anisotropy is biaxial with easy 
axes aligned with the <110> pc lattice direction o 9 ' 10 . 
(We use subscript c, pc, o and t for cubic, pseu- 
docubic, orthorhombic and tetragonal crystal struc- 
tures, respectively.) Next to the magnetocrystalline 
anisotropy a uniaxial anisotropy is present as well, 
which is stepedge induce d 13 ! 16 . Here we investigate the 
case of isotropic compressive strain, which can be real- 
ized with LSMO thin films on the cubic (LaA103)o.3- 
(Sr 2 AlTaO 6 ) .7 (LSAT) (001) c substrate. LSMO thin 
films under compressive strain adopt an orthorhombic 



crystal structur e) 17 ! 18 , which is characterized by the pres- 
ence of oxygen octahedra rotations around all three pseu- 
docubic crystal axes. As the magnetic coupling depends 
on the Mn-O-Mn bond angl o 19 ' 20 , it is an interesting 
question whether the magnetic properties are anisotropic 
in the different orthorhombic directions. Note that for 
another case, orthorhombic LSMO grown on NGO (110) o 
the difference in lattice mismatch between the two in- 
plane directions determines the anisotropy^, so this sys- 
tem is not suitable to study the effect of the orthorhom- 
bicity on the magnetic properties. For LSMO films grown 
on NGO (110)o the [110] o lattice direction is subjected 
to less compressive strain than the [001] o lattice direction 
and is therefore the easy axis due to the strain anisotropy. 
For LSMO films grown on LSAT the lattice mismatch is 
equal and the anisotropy is due to the intrinsic anisotropy 
of the orthorhombic crystal structure between the [110] o 
and [001] o lattice directions. 



Here, we show that LSMO thin films can be grown 
coherently and untwinned on LSAT substrates and that 
the orthorhombicity induces anisotropic magnetic prop- 
erties. Next to a biaxial component of the magnetic 
anisotropy, we observed a uniaxial component to the 
anisotropy which is aligned with the principal crystal di- 
rections and became more pronounced for increasing film 
thickness. We found no correlation between the uniax- 
ial anisotropy and the stepedge direction. We obtained 
twinned samples, by growth on surfaces with reduced 
crystalline quality, for which the uniaxial anisotropy 
was reduced. Therefore we conclude that the uniaxial 
anisotropy is caused by the orthorhombic crystal struc- 
ture. 
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FIG. 1. (Color online) Surface analysis of the LSAT substrate 
by atomic force microscopy, a) after annealing at 1050° C for 
12 hours, b) after annealing at 950° C for 1 hour, The images 
are 5 by 5 fim and the color scale is 2 nm. The insets show a 
close-up of the roughness of the terraces. 



II. SAMPLES AND SUBSTRATE 
PREPARATION 

The as-received LSAT substrates were cleaned with 
acetone and ethanol before they were subjected to an an- 
neal treatment. Two anneal treatments were used to ob- 
tain respectively surfaces with smooth terraces and sur- 
faces with sub unit cell roughness on the terraces. The 
first treatment consisted of an annealing step at 1050°C 
for 12 hour in 1 bar of O2 gas pressure. For the second 
treatment both the anneal time and temperature were 
decreased to 1 hours and 950° C respectively. The sur- 
faces were characterized with atomic force microscopy 
(AFM) . Typical results are shown in figure Q] For the 
substrates subjected to the first anneal treatment a step 
and terrace structure with 4 A (a single unit cell) step 
height was observed. The stepedges were not straight 
but meandering and 4A deep holes are observed near the 
stepedges. Note that the miscut of these substrates is 
very small, approximately 0.02°, giving a terrace width 
of more than 1 /im. Between the stepedges areas with 
atomically smooth morphology were observed. The sub- 
strates subjected to the second treatment show terraces 
with reduced crystalline quality, but still single unit cell 
step heights. 

LSMO thin films were grown on the LSAT (001) sub- 
strates by pulsed laser deposition (PLD) from a stoichio- 
metric target in an oxygen background pressure of 0.35 
mbar with a laser fluence of 3 J/cm 2 and at a substrate 
temperature of 750°C. After LSMO deposition, the films 
were cooled to room temperature at a rate of 10°C/min 
in a 1 bar pure oxygen atmosphere. The growth settings 
were previously optimized and were identical to the ones 
used for LSMO films on other substrate s 14 ' 21 . 

In this paper four samples are described, see table |U 
Sample U12 and U40 were grown on substrates with a 
smooth surface and have a thickness of 12 and 40 nm 
respectively. Samples T29 and T50 were grown on sub- 
strates with terraces with reduced crystalline quality and 
are respectively 29 and 50 nm thick. (The sample labels 
consist of either the letter T or U for twinned/untwinned 



and a number which indicates the sample thickness.) The 
sample thicknesses were measured with x-ray reflectiv- 
ity measurements. AFM measurements (not shown) re- 
vealed surfaces of the thin film where the morphology of 
the substrate was still visible. The Curie temperature of 
the films was larger than 350 K (350 K was the measure- 
ment limit of the vibrating sample magnetometer) and 
did not depend on film thickness and the twinning, as 
discussed in the next section, of the films. 



III. STRUCTURAL CHARACTERIZATION 

The top panel of figure [2] shows reciprocal space maps 
of LSMO and LSAT around the (204) c , (024) c , (204) c 
and (024) c LSAT reflections. These results were ob- 
tained from sample U40. The LSMO has a slightly dis- 
torted orthorhombic (monoclinic) unit cell with (110) o 
out-of-plane orientation and (110) o and (001) o in-plane 
orientations. The orthorhombicity can be deduced from 
figures [2k and which show a difference in lattice 
spacing for the (260) o and (620) o LSMO reflections, 
which represent a dissimilarity between the LSMO a 
and b lattice parameters. The lattice parameters are as 
follows: a=5.47±0.0lA, b=5.51±0.0lA, c=7.74±0.0lA, 
a=90±0.1°, /3=90±0.1° and 7 =89.6±0.1°. 

Next to the LSMO (444) and (444) reflections satel- 
lites are observed. In a previous paper— we have dis- 
cussed the orthorhombic crystal structure and the pres- 
ence of satellite peaks in the [001] direction for LSMO 
grown on NGO (110) o substrates. The satellites result 
from periodic lattice modulations which partially relieve 
the applied strainer—. As the LSMO crystal structure 
can easily relieve strain in the [110] o lattice direction with 
a change in the 7 angle, the lattice modulations are only 
present in the [001] o direction. We conclude that LSMO 
films on LSAT behave similarly as LSMO on NGO. 

The panels e-h) of figure [5] show the same reciprocal 
space maps, but obtained for sample T50. Dissimilar 
lattice spacings are observed in both in-plane directions 
and the satellite peaks are visible in all reciprocal space 
maps. Zhou et al. observed similar satellite peaks around 
reflections of LSMO on LSAT and attributed the satel- 
lites to an in-plane superlattice with alternating [110] o 
and [Tl0] o orientations 2 ^.. This cannot be the case for 
our samples as the [Tl0] o in-plane orientation would re- 
sult in reciprocal space maps with both the (620) o and 
the (260) o peak visible in the same plot. In e.g. figure [2f 
this is clearly not the case. Together with the presence of 
the satellites in figure [2Jd and[2]i only in the [001] o direc- 
tion, we conclude that sample D is twinned in domains 
with different [110] o and [001] o orientations. 

For the two thinner samples we observed the same be- 
havior as the thicker samples with equal substrate treat- 
ment (not shown). Sample U12 was untwinned as con- 
cluded from the positions of the Bragg reflections, but 
here no satellites could be observed. Sample T29 was 
twinned with satellites in both directions in reciprocal 
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Sample 


Thickness 


Substrate 


Crystal 


Satellites 


0easy ( ) 


k u /ki 






M 


surface 










(J/m 3 ) 


U12 


12 


smooth 


untwinned 


not observed 


40±1 


0.18 ±0.05 


110±30 


U40 


40 


smooth 


untwinned 


along [001] o 


12±1 


0.91 ±0.03 


540±5 


T29 


29 


rough 


twinned 


both directions 


42±1 


0.1 ±0.05 


60±30 


T50 


50 


rough 


twinned 


both directions 


31±1 


0.48 ±0.03 


290±25 



TABLE I. The ratio between the anisotropy constants for the various samples. The angle of the easy axis is obtained from the 
fitting procedure and the ratio between the anisotropy constants is calculated with equation [2] 



space. 

To explain the absence of twinning in the samples with 
smooth surfaces we compare our results to the more 
widely studied SrRu0 3 (SRO) thin films on STO sub- 
strates. SRO is orthorhombic and single domain films can 
be realized by growth on smooth vicinal substrates with 
stepedges approximately aligned with the main crystal 
axis. In that case the [001] o lattice directions aligns 
with the stepedge direction^. This has been explained 
in three different ways. Gan et al. suggested that single 
domain growth is related to the observed stepflow growth 
mode2&. Maria et al. suggested that stepedge strain is 
the dominant mechanism as the films are not orthorhom- 
bic at deposition temperature^. Finally Vailionis et al. 
suggested that the films are tetragonal at deposition tem- 
perature and that the [001]t lattice direction aligns pref- 
erentially with the stepedges^. In contrast to SRO the 
LSMO films have a preferred orientation with the [001] o 
lattice direction aligned perpendicular to the stepedges 
and does not grow in stepflow mode. Therefore the ex- 
planation by Maria et al. is the most likely candidate 
to explain the single domain growth. During cooldown 
orthorhombic domains nucleate at the stepedges and the 
stepedge strain favors octahedra buckling in one direction 
and the domains continue to grow across the terraces re- 
sulting in a single domain film. For the samples with 
a rough surface the orthorhombic domains can nucleate 
at defects at the substrate surface and no preferential 
orientation exists, resulting in twinned films. 

In summary, growth on a relatively smooth surface re- 
sults in untwinned LSMO films, while growth on terraces 
with reduced crystalline quality results in twinned LSMO 
films. We used the difference in magnetic properties be- 
tween twinned and untwinned films to identify the contri- 
bution from the orthorhombicity as it should be reduced 
for the twinned samples. 



IV. MAGNETIC CHARACTERIZATION 

The samples were characterized with vibrating sam- 
ple magnetometer (VSM, Model 10 VSM by Microsense) 
measurements at room temperature. The in-plane angle 
of the applied field was varied to determine the mag- 
netic anisotropy. For all field angles a full magnetization 
loop was measured and the remanent magnetization was 
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FIG. 2. [Color online] Top panel: Reciprocal space maps 
around the a) (204) c , b) (024) c , c) (204) c and d) (024) c LSAT 
reflections of the 40 nm thick sample grown on a smooth sub- 
strate (sample U40). In a) and c) the dissimilar spacing of 
the (260) o and (620) o LSMO reflections is clearly observed 
while in b) and d) satellites are present next to the (444) 
and (444) LSMO reflections. Middle panel: Reciprocal space 
maps around the e) (204) c , f) (024) c , g) (204) c and h) (024) c 
LSAT reflections of the 50 nm thick sample grown on a rough 
substrate (sample T50). The satellites are present in all maps 
and both e) and h) show intensity at the position for the 
(260) o LSMO reflection. Therefore this sample is twinned. 
Bottom panel: Schematic of the real space crystal structure 
of an untwinned sample viewed along i) [001] o and j) [110] o - 
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obtained from the loop. Figure [3^ shows the magneti- 
zation loops of sample U12 with the field aligned with 
three high symmetry directions. In figure we plotted 
the dependence of the remanent magnetization on the in- 
plane field angle. The largest remanence was found for an 
applied field at approximately 40 degrees with respect to 
the [001] o in-plane lattice direction. A predominant biax- 
ial behavior is observed with easy axes aligned with the 
[110]pc and symmetry related crystal directions. Next 
to this biaxial anisotropy a small uniaxial anisotropy is 
present as well which can be seen in the difference in re- 
manent magnetization at and 90 degrees and the shift 
of the easy axes to ±40°. This uniaxial contribution be- 
comes more pronounced in the thicker film (sample U40) 
shown in figure [3h and[3ji. The remanent magnetization 
of sample B at degrees approaches the easy axes value 
and the remanent magnetization at 90 degrees is much 
smaller, only 20 percent of the easy axes value. Due to 
the combination of uniaxial and biaxial anisotropy the 
easy axes are shifted to ±15°. The hard axis magnetiza- 
tion loop does not show a switch of the magnetization, 
but an almost linear dependence of the magnetization on 
the applied field which is characteristic for a hard axis 
loop in a sample with a uniaxial anisotropy. 

A similar thickness dependence of the anisotropy was 
also found for the samples T29 and T50, which are 
twinned. The results are plotted in figure E^-Oi- A biax- 
ial and a uniaxial contribution are present and the uni- 
axial contribution is more pronounced in the thicker film. 
Sample T50 shows easy axes which are shifted to ±30° 
and the hard axis remanence is 50 percent of the easy 
axis value. Comparing samples U12 and U40 with sam- 
ples T29 and T50 we find that the uniaxial contribution 
is more pronounced in the samples U12 and U40. 



V. DISCUSSION 

In order to quantify the biaxial and uniaxial contribu- 
tion to the anisotropy we start with a general anisotropy 
energy equation which contains both a biaxial and a uni- 
axial contributior 
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E a /V = k n sin 2 



+ |sin 2 (2(^- 
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in which E a is the anisotropy energy, V the volume of 
the sample, fc u (fci) the uniaxial (biaxial) anisotropy con- 
stant, <fi the angle of the magnetization, 0i=O° (</>2=45°) 
the angle of the easy axis of the uniaxial (biaxial) 
anisotropy. The easy axes are found by minimizing the 
energy with respect to <f>. This results in: 
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k-a/ki for fc u < k\, 
for fc u > fci, 



(2) 



from which the easy axes, 4> eaS y> can be obtained. The 
measured remanence versus field angle dependencies are 
the projections of the magnetization, which at remanence 
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FIG. 3. [Color online] a) Magnetization loops with the field 
aligned with 3 high symmetry directions and b) remanence 
versus field angle dependence obtained from sample U12. The 
latter graph shows the data from the measurements (black 
dots) and the result from the fit procedure described in the 
text (blue line), c) and d) The corresponding graphs for sam- 
ple U40. e) and f) Sample T29. g) and h) Sample T50. The 
field angle is defined with respect to the crystal structure, 0° 
(90°) corresponds to the [001] o ([110] o ) lattice direction. 



is aligned with the easy axis, onto the measurement di- 
rection: 



M rem (0) = M cos((9 



(3) 



where M rom is the remanent magnetization, 9 is the field 
angle, Mq is the remanence in the easy direction and (/> eaS y 
is the closest easy axis. Equation [3] has been used to ob- 
tain the fits in figure[3] The measured data, and therefore 
the anisotropy, is well described by equation Q30 This 
allows us to extract the ratio between the uniaxial and 
biaxial anisotropy energies k u /ki. The results are pre- 
sented in table [I] For all samples the uniaxial anisotropy 
energy is found to be smaller than the biaxial anisotropy 
energy. Between the different samples the ratio k u /ki 
changes by a factor of 10. 

The anisotropy field H an can be obtained from the 
slope dM j dH at H =0 in the hard axis magnetization 
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loop of a material with uniaxial anisotropy by the rela- 
tion 



M a , 



dM/dH ' 



(4) 



Here, M sat is the saturation magnetization, M is the 
magnetization and H is the applied field. The anisotropy 
energy is given by 



2fc u = # an /xoM s: 



(5) 



in which /io is the permeability of free space. We obtained 
a value for k u of 540±5 J/m 3 for the almost uniaxial sam- 
ple U40. This would imply a biaxial anisotropy constant 
of 600±10 J/m 3 which corresponds well with earlier ob- 
tained results2i2i. Therefore one can assume that the k\ 
value is the same for all films and independent of thick- 
ness. The uniaxial anisotropy constants are calculated 
from the k u /ki ratios and presented in table U as well. 
The uniaxial anisotropy is very weak for these samples, 
only in the range 50-600 J/m 3 . 

Next we turn to the origin of both contributions to 
the anisotropy. The biaxial contribution with easy axes 
aligned with the <110> pc lattice directions corresponds 
well with earlier results of magnetic anisotropy of LSMO 
on STO (001)c substrates*^. This represents the intrin- 
sic magnetocrystalline anisotropy of LSMO strained to 
in-plane cubic symmetry. For the uniaxial contribution 
different explanations exist. It is well known that a weak 
uniaxial anisotropy in LSMO can be the result of ste- 
pedge induced anisotrop y 13 ' 16 . However, the thickness 
dependence of the uniaxial anisotropy is at odds with an 
interpretation in terms of stepedge induced anisotropy. 
The stepedge induced anisotropy compared to the biax- 
ial anisotropy should scale with the ratio of the volume 
of the surface layers containing the stepedges (in practice 
the monolayers at the interfaces) to the volume of the 
film. Therefore the stepedge induced anisotropy should 
be most pronounced for the thinnest films as the miscut 
of the samples was comparable. Here the opposite is ob- 
served. Also the uniaxial easy axis was not aligned with 
the stepedge directions in most of the samples, sometimes 
even 90 degrees perpendicular to the stepedges. This 
rules out the contribution of stepedge induced anisotropy. 

The uniaxial easy axis is aligned with the [001] lattice 
direction while the hard axis is aligned with the [110] o lat- 
tice direction. This, together with the reduced uniaxial 
anisotropy for the twinned samples, relates the observed 
anisotropy to the orthorhombicity of the LSMO films. It 
is unclear what the origin of the magnetic anisotropy in 
the orthorhombic crystal structure is. In general mag- 
netic anisotropy in manganites is explained by a global 
strain field which relates the easy axis to the maximum 
strain direction 3 ^. This cannot be applied for these sam- 
ples, as the LSMO has equal strain in the [001] o and the 
[110] o lattice directions due to the cubic symmetry of the 
substrate. Note that for the case of orthorhombic LSMO 
on NGO (110) o the uniaxial easy axis is aligned with the 



[110] o lattice direction and is due to strain anisotropy 14 , 
in contrast to the observed anisotropy of LSMO on LSAT. 

The difference between the two lattice directions in the 
orthorhombic symmetry is due to the different oxygen 
octahedra rotations. We therefore suspect that the mi- 
croscopic origin of the magnetocrystalline anisotropy in 
LSMO is effected by the oxygen octahedra rotations. A 
plausible scenario would be that the difference in octahe- 
dra rotations along the [110] o and [001] o directions results 
in a different Mn-O-Mn bond angle (< Mn-O-Mn >)in 
these two directions. For manganites it is well known 
that the bandwidth, due to the double exchange mech- 
anism, depends on the Mn-O-Mn bond angle according 



W = 



d 3.5 

u MnO 



(6) 



in which W is the bandwidth, u = 0.5(7r— < 
Mn-O-Mn > is proportional to the amount of octahedra 
rotation and dMnO is the Mn-0 distance. This suggests 
that for LSMO on LSAT the Mn-O-Mn bond angle is 
larger in the [001] o direction than in the [110] o direction, 
resulting in an increase in the bandwidth and therefore 
the easy axis is aligned with this direction. This explana- 
tion is also consistent with the increase of k u /k\ with film 
thickness. Due to structural reconstructions at the inter- 
face and the surface 3 - 1 -, the oxygen octahedra rotations in 
the surface and interfacial regions deviate from those of 
the bulk of the film. The biaxial anisotropy is intrinsic 
to MnOg octahedra and not sensitive to the structural 
reconstructions. In a follow up work we will investigate 
the crystal structure, and especially the octahedra rota- 
tions, in more detail, in order to resolve the microscopic 
origin of the magnetocrystalline anisotropy in LSMO. 

An alternative interpretation of the magnetic data is 
that the anisotropy is somehow related to the lattice 
modulations observed with the satellites in the XRD 
measurements. As the lattice modulations only occur 
in the orthorhombic [001] direction, it is not possible 
to discriminate between anisotropy induced by the or- 
thorhombic crystal structure and anisotropy induced by 
the lattice modulations in our experiment. Neverthe- 
less, we expect that the lattice modulations result in mi- 
crotwins which have reduced magnetic coupling at the 
microtwin boundaries. The shape anisotropy of each mi- 
crotwin would then be aligned with the [110] o direction, 
contrary to the observation of a magnetic [001] easy axis. 

Although one would suspect that the uniaxial 
anisotropy would disappear for the twinned samples, this 
is not the case. We assume that this is due to the domi- 
nant presence of grains with one orientation. 



VI. CONCLUSION 

LSMO films with an orthorhombic crystal structure 
can be grown coherently and untwinned on cubic LSAT 
substrates. The magnetic anisotropy of the films is 
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described by a combination of biaxial anisotropy with 
easy axes along the <110> pc directions and a uniaxial 
anisotropy with easy axis along the [001] o direction. For 
thicker films the uniaxial anisotropy becomes more pro- 
nounced. The uniaxial part of the anisotropy is induced 
by the orthorhombic symmetry of the LSMO. We expect 
these findings to lead to a better understanding of the 



microscopic origin of the magnetocrystalline anisotropy 
in LSMO. 
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